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Abstract. We present XMM-Newton MOS imaging and PN timing data of the millisecond pulsar PSR J02 18+4232. We 
confirm the previously detected pulsations of PSR J02 18+4232 and we show that the folded lightcurve is dependent on energy. 
We present the broad band (0.2-10.0 keV) spectrum of this millisecond pulsar, as well as the spectra of: the pulsed emission, 
the individual pulses, the interpulse and the non-pulsed region and we compare our results with previous data from Rosat, 
BeppoSAX and Chandra. We discuss the results of the spectral fitting in the context of recent pulsar models. 
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1. Introduction 

Pulsars show two different types of X-ray pulsations: soft 
pulsed emission with blackbody spectra where the emission is 
from the neutron star surface, such as cooling emission from 
the whole surface as a relic of the neutron star's birth and/or 
the emission from a heated region, most likely the polar cap; 
or pulsed emission with narrow pulses and hard power law 
spectra, where it is tho ught that the emission arises from the 
magnetosphere (see e.g . lSaito et al. Hl997UKuiper et al. fl 998: 
iTakahashi et al. 1 1200 lL for further discussion). Therefore de- 
tecting not only a pulsation, but also determining the form 
of the pulse profile, could help to identify the underlying 
emission mechanism and detection of the polar cap ther- 
mal radiation would allow us to discriminate between var- 
ious models of p ulsars and study the properties of neutron 
star surface layers JZavlin et al. 1120021 IPavlov & Zavlinll997l 
IZavlin & Pavlov I fl998l) . Many millisecond pulsa rs (MSPs) 
such as PSR B 182 1-24 (e.g. ISaito et al. I M 997l) and PSR 
B1937+21 (e.g. Takah ashi et al. I l200lh . like many ordinary 
pulsars, show X-ray pulsations and have been observed to 
emit X-ray spectra that are well fitted by a hard power law. 
In contrast, it has been proposed that the millisecond pulsar 
PSR J04 37-4715 emits thermal radiation fr om hot polar caps 
jZavlin & Pavlov 1 1998tlZavlin et al. l20ol . 

Th e radio source J0218+4232 has been known for many 
years JOwarakanath A Shankarlfl99rl lHales et al. lfl99l al- 
though it was first c onfirmed to be a millisecond pulsar by 
Navar ro et al. I ll 19951) . using radio observations made with the 
Lovell telescope in 1993. A pulse period of 2.3 ms was deter- 
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mined. iNavarro et all dl995h also showed that this luminous 
pulsar (L400 > 2700 mJy kpc 2 lNavarro et al. 1 19951 where L400 
is the luminosity at 400 MHz) is in orbit with a low mass white 
dwarf ( O.2M ), with an orbital period of about two days. From 
the dispersion measure a lower limit on the distance of 5.7 kpc 
was derived. 

Verbu nt et al. ] Jl996) showed evidence for PSR 
J0218+4232 at high energies, both in X-rays using ROSAT 
and y-rays, using EGRET, the High-Energy Gamma-Ray 
Telescope aboard the Compton Gamma-Ray Observatory 
(CGRO). They also present ed some indication of pulsation in 
both energy domains. Kuip er et al. I d2000l) provided stronger 
evidence for the detection of pulsed y-ra y emission fro m PSR 
J0218+4232, also using EGRET data. Kuiper et al. I (fl998) 
presented further evidence for X-ray pulsation at the radio 
pulse period using 98 ks of ROSAT data. They also likened 
the sharp pulses to those of PSR B 182 1-24 and the Crab 
pulsar, which indicates a magnetospheric origin of the pulsed 
X-ray emission. They also estimated that almost two thirds of 
the soft X-ray (0.1-2 .4 keV) emission is non-pulsed, where 
INavarro et al. I i 19951) also stated that a large fraction of the 
radio emission is not pulsed. This lead them to the conclusion 
that PSR J02 1 8+4232 is a n aligne d rotator. This was further 
supported by Stairs et al. I (E?99), using polarimetric radio 
observations. 

BeppoSAX o bservations of PSR J02 18+4232 
dMineo et al. I l2000l) provided detailed information on the 
pulsar's temporal and spectral emission prope rties (1.6-10 
keV). They showed that pulse 1 (lying at 0=0.8, in M ineo et al. I 
2000) is much stronger than pulse 2 (lying at 0=0.3) in the 
1 .6-4 keV energy band, but that pulse 2 is much stronger than 
pulse 1 in the 4-10 keV energy band. Also at the higher energy 
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no non-pulsed component was apparent above the background 
level. They determined a very hard power law fit, photon 
index = 0.61+0.32, to the pulsed part of the spec trum, with 
a fixe d column density of 5xl0 20 cirT 2 (see IVerbunt et al. I 
Il996l) . which gave a flux of 4.1xl0~ 13 ergs cm" 2 s" 1 (2-10 
keV). The measured spectral indices of pulse 1 and pulse 2 
were 0.84+0.35 and 0.42+0.36. The total spectrum was well 
fitted by a power law, with a photon index of 0.94+0.22 and 
they found a flux of 4.3xl0~ 13 ergs cirT 2 s" 1 (2-10 keV). 
Kuio er et al. I fco02a) showed that the non-pulsed spectrum 
of PSR J0218+4232 has a softer spectrum than the pulsed 
emission, usin g 73 ks of Ch a ndra data . 

Recently, Kuiperetal. ] J2002albh presented multiwave- 
length (radio, X- and y-ray) pulse profiles of PSR J0218+4232. 
They found that the two X- and y-ray pulses are aligned and 
that the X-ray pulses are also aligned with two of the three ra- 
dio pulses. This is a very interesting result as multiwavelength 
emission can help to discriminate between the different theoret- 
ical models proposed to describe the nat ure of m illisecon d pul- 
sars, such as the polar cap models (see lHarding & Muslimovl 

2002 1 and referen ces therein) and the outer gap models (see 

~ "™— it -i ° r 

Cheng et al. 2000, and references therein), and thus trace the 

origin of the emission. For example, in polar cap scenarios, y- 
ray pulses should be accompa nied by X-ray pulses at the same 
phase and with similar shapes (Rudak & Dvks 1999). 

XMM-Newton combines the spectral capabilities of 
ROSAT (0.1-2.5 keV) and the BeppoSAX MECS (2-10 keV), 
which, when coupled with the XMM-Newton's sensitivity, 
gives us unparalleled spectral information from low to high en- 
ergies. We have thus used XMM-Newton to observe the mil- 
lisecond pulsar PSR J0218+4232 to show the evolution of the 
pulsed profile from low energies (0.4-1.6 keV) to high energies 
(up to 12 keV). We present spectra of the pulsar as a whole 
and both the pulsed and non-pulsed components and we draw 
some conclusions about the nature of the emission from PSR 
J0218+4232. 

2. Observations and data reduction 

PSR J0218+4232 was observed by XMM-Newton on 2002 
February 11-12. The observations spanned 37.2 ks (MOS cam- 
eras) and 36 ks (PN camera), but a soft proton flare affected 
17 ks of the MOS exposure and 16.5 ks of the PN expo- 
sure. The MOS data were reduced using Version 5.3.3 of 
the XMM-Newton SAS (Science Analysis Software). However, 
for the PN data we took advantage of the development track 
version of the SAS. Improvements have been made to the 
OAL (ODF (Observation Data File) access layer) task (version 
3.106) to correct for spurious and wrong time values, prema- 
ture increments, random jumps and blocks of frames stemming 
from different quadrants in the timing data in the PN auxil- 
iary file, as well as correcting properly for the onboard delays 
(Kirs ch etal. 1 12003). Using this version, our timing solution 
improved (see Sect. [3. 

We employed th e MOS cameras in th e full frame mode, us- 
ing a thin filter (see iTurner et al. 1 l200ll) . The MOS data were 
reduced using 'emchain' with 'embadpixfind' to detect the bad 
pixels. The event lists were filtered, so that 0-12 of the pre- 



defined patterns (single, double, triple, and quadruple pixel 
events) were retained and the high background periods were 
identified by defining a count rate threshold above the low 
background rate and the periods of higher background counts 
were then flagged in the event list. We also filtered in energy. 
We used the energy range 0.2-10.0 keV, as recommended in 
the document 'EPIC Status of Calibration and Data Analysis' 
( Kirs ch et al. 1120021) . The event lists from the two MOS cam- 
eras were merged, to increase the signal-to-noise. 

The PN camera was also used with a thin fil ter, but in tim- 
ing m ode which has a timing resolution of 30jus JStruder et al. I 
2001). The PN data were reduced using the 'epchain' of the 
SAS. Again the event lists were filtered, so that 0-4 of the 
predefined patterns (single and double events) were retained, 
as these have the best energy calibration and we filtered in 
energy. T he document ' EPIC Status of Calibration and Data 
Analysis' ( Ki rsch et al. 1120021) recommends use of PN timing 
data above 0.5 keV, to avoid increased noise. We used, in gen- 
eral, the data between 0.6-12.0 keV as this had the best signal- 
to-noise, although we have also tried to exploit the data in the 
0.4-0.6 keV band (see Sect. [3). The on-board event times ex- 
pressed in the local satellite frame were subsequently converted 
to Barycentric Dynamical Time, using the SAS task 'barycen' 
and the coordinates derived f rom 6.5 years o f radio timing mea- 
surements, as given in Kuiper et al. (20021). 

3. Timing analysis 

In the PN timing mode, all the two-dimensional spatial infor- 
mation is collapsed into a single dimension (x-direction). We 
extracted the pulsar data using a rectangular region centered on 
the pulsar. 

We corrected the data for the orbital motion of the pul- 
sar and the d a ta were folded on the radio ephemeris given in 
IKuiper et al. I (|2002a). We verified that the corrections made 
to the timing data with the development version of the SAS 
are indeed correct. We tested frequencies at and around the 
expected frequency. We used a x 2 test t0 determine the fre- 
quency at which the pulse profile was the strongest. We found 
the largest peak in the xl versus change in frequency from the 
expected frequency at Av = 2 xl0~ 6 s~', see Fig. ^ This is 
well inside the resolution of this dataset (~1/T ^ J ), which is 
3 xl0" 5 s _1 , thus we can conclude that the data reduction and 
analysis made to the dataset are reliable. The folded lightcurve 
(0.6-12 keV), data folded on the radio frequency, counts ver- 
sus phase, is shown in Fig. |2] This lightcurve is ver y simi- 
lar to the lightcurve presented bv lKuiper et ail (|2£)02a), where 
both pulses can be found at approximately the same phase 
(within the ti ming uncertainties of Chandra and X MM-Newton 
IKuiper et al. ll2002aHTennant et al. lEoOlllKirsch et al. 112003b . 
Fitting the two pulses wit h two Lorentzians a nd a back- 
ground (in a similar way to Ku iper et al. Il2002al) . we find a 
pulse separation of 0=0.491+0.022 (9 0% confidenc e limit) , 
consistent with pr evio us results e.g. IKuiper et al. I (^998), 
iMineo et al. I feOOOl) and lKuiper etaH l l2002al) . The FWHM of 
the pulses are <S0i =0.1 12+0.038 for the pulse centred at phase 
01=0.242+0.008 and 502=0.121 +0.056 for the pulse centred at 
02=0.733+0.014. These values are also consistent with the val- 
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Fig. 1. xl versus change in frequency from the expected pulsa- 
tion frequency (shown as the soli d vertical l i ne at Ay = 0.0), cal- 
culated with the ephemeris o flKuiper et al. Id2002al) . See Sect.[3] 
for details. 

ues presented in Kuip er etal. 1 12002a). Fitting the pulses with 
two Gaussians or fitting the pulses separately gives a similar 
result. 

To calculate the DC (non-pulsed) level of our lightcurve 
we us ed the bootstrap method outlined by Swan epoel et al. I 
Jl996l) . The DC level is shown in Fig.|2]by the solid line and the 
two dashed-dotted lines show the + l<x errors. The background 
level (also shown in Fig. [2 dashed line) was determined by ex- 
tracting the data from a region in the same observation as the 
pulsar, centered 12 pixels (~50") in the x-direction from the 
pulsar, where there were no other X-ray sources close to this 
position or in the full length of the y-direction. We verified that 
there were indeed no sources in this region, using the MOS 
imaging data (see Section|4}. The extraction region chosen for 
the background was the same shape and size as the extraction 
region for the pulsar. We see evidence for both a DC inter- 
val between = 0.45-0.58 (see Fig. |2} and an interpuls e be- 
tween <b = 0.88-1.1 1, in similar places to those of lKuiper et al. I 
(2002a), where the interpulse region has 10.5+3.3 counts/bin 
more than the non-pulsed region. Thus we conclude that our 
X -ray pulse profile can be compared directly to the lightcurve 
of lKuiper et al. I J2002al) . We therefore refer to the pulse situ- 
ated at 0i=O.24 as pulse 1 (phase = 0.1 1-0.45) a nd <fo=0.73 as 
pulse 2 (phase = 0.58-0.88) in the same way as IKuiper et al. I 
(2002a). Subtracting the background from the total count dis- 
tribution measured in the strip centred on the pulsar, we found a 
pulsed percentage of 69+6% in the 0.6-12 .0 keV band, similar 
to the value determined by Kui per et al. I (l2002al) . who found 
64+6% with Chandra in the 0.08-10.0 keV band. We find that 
pulse 1 contributes 65 ±7% of the pulsed percentage and pulse 
2 contributes 32+8%, where the ratio between the pulse 1 and 
the pulse 2 counts is 1.60+0.09 in the 0.6-12.0 keV band. 

Given the large number of counts in our observations, we 
were also able to analyse the variations in the lightcurves as a 
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Fig. 2. Lightcurve folded on the ephemeris of IKuiper et al. I 
(2002a) and binned into 30 bins, each of 77/zs. Two cycles 
are shown for clarity. Pulses 1 (PI) and 2 (P2) are indicated. 
Between these two pulses is the section defined as the non- 
pulsed region (NP). Between pulses 2 and 1 is the section de- 
fined as the interpulse (IP). A typical ±l<x error bar is shown. 
The dashed line shows the background level, where the error 
bar represents the + l<x error. The solid line shows the DC level 
and the dashed-dotted lines represent the + l<x error. 



function of energy. We chose the energy bands: 0.4-1.6 keV; 
1.6-4.0 keV; and 4 .0-10.0 keV, so that we could compare our 
results to those of iMineo et al. I d2000l) . see Fig. [3] The data 
were folded as before and binned into 1 8 bins, where two cycles 
are presented for clarity. 

The first thing that we notice from looking at Figgis that 
the lightcurves are remarkably different in the 3 energy bands, 
indicating that the spectra of the pulses are indeed quite dissim- 
ilar. Pulse 1 is always the strongest of the two pulses, but the 
ratio between the number of counts in the two pulses in the dif- 
ferent energy bands varies quite dramatically. This can be seen 
quantitatively in Tabled The softest and intermediat e band s 
are consi stent with th e ratio s found bv IMineo et al. I (|2000). 
However Min eo et al. I ([2000) found that pulse 2 was stronger 
than pulse 1 in the hardest band (0.4-10.0 keV), but they use 
a pulse width of 0=0.37 for pulse 2, whereas we use a pulse 
width of 0=0.30, which could account for the discrepancy. We 
also see evidence for the interpulse in the softest energy b and. 

As the ROSAT and Be ppoSAX phasograms ( Kuipe r^t al. I 

Il998t IMineo et al. I l2000l) are not in absolute phase, it is 
important to deter mine which of the pulses correspond to 
our pulse 1 and 2. IKuiper et al. present the tentative 

alignment after cross-correlation of the two X-ray lightcurves 
(ROSAT HRI, 0.1-2.5 keV and BeppoSAX MECS 1.6-10.0 
keV), where pulse 1 (the largest pulse in both the ROSAT and 
the BeppoSAX lightcurves) app ears at approx imately 0=0.2, 
as in our lightcurve and that of lKuiper et al. I J2002al) . which 
is in absolute phase. Thus we can compare pulses 1 and 2 in 
IMineo et al. I (yOOO) with pulses 1 and 2 in this work. 
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Table 1. Table showing the pulsed percentage and the ratio of 
pulse 1 to pulse 2 for the PSR J02 18+4232 in different energy 
bands. 



Fig. 3. Lightcurves folded on the pulse period and binned into 
18 bins, (a) 0.4-1.6 keV (b) 1.6-4.0 keV (c) 4.0-10.0 keV. Two 
cycles are shown for clarity. The dashed line indicates the back- 
ground level, the dashed-dotted line indicates the DC level and 
an error bar represents the typical +lcr error. 



Comparing our lightcurves with those presented in 
iMineo et al. I d2000l) . we find a remarkable visual difference. 
In the intermediate energy band (1.6-4.0 keV), they find only 
one clear peak - pulse 1 . However, we find that we have two 
pulses of similar strength. In the hardest band, we see that 
pulse 1 is strongest, in c ontrary to the BeppoSAX observa- 
tions dMineo et al. Il200ol) . The discrepancies have yet to be 
explained, but they may be related, to some extent, to differ- 
ences in the energy response of the different instruments. In the 
softest band (0.4-1.6keV), we see that pulse 1 is strongest, sim- 
ilar to ROSAT(0.1-2 4keV) JKuiper et al. ll998HMineo et al. I 
l2000l) . lKuiDer et al. I (E)02b) also find evidence that pulse 2 is 
stronger at higher energies than pulse 1, using RXTE data in 
the range 2-16 keV, as well as some evidence that pulse 2 is 
stronger than pulse 1 in the 2-8 keV range. Examining the 8-12 
keV range, we find that pulse 2 is 1.5 times stronger than pulse 
1. From this we have seen that the two pulses are spectrally 
quite dissimilar. With this in mind we investigated the pulse' 
spectra, as well as other regions of the lightcurve (see Sect.|5j. 



4. The MOS cameras' field of view 

The sources in the MOS field of view were first detected us- 
ing the SAS EPIC source detection task 'eboxdetect', which 
employs a 'local' source detection algorithm. A box of 5 x 5 
pixels was used to detect point sources and then the same box 
was used on the background. To detect extended sources, two 
iterations were made with a box of 10 x 10 pixels and then 
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0.6-12.0 


69±6 


1.60+0.1 


0.4-1.6 


59+7 


2.1+0.1 


1.6-4.0 


68±11 


1.1+0.1 


4.0-10.0 


43±12 


2.5+0.2 
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Fig. 4. MOS (20.5 ks) image of the pulsar field of view. The 
pulsar can be seen as the brightest source in the centre of the 
field of view. 

20 x 20 pixels, however only point sources were detected. A 
detection likelihood could then be calculated. 

42 sources have been detected using the SAS task 'emlde- 
tect', using a maximum detection likelihood of 18 (approxi- 
mately 4<x likelihood of detection) and ignoring those sources 
not found on both cameras (unless they lay outside the field of 
view or were in a chip gap). We calculate an unabsorbed flux 
limit of 1.0 xlO~ I4 ergs cirrV 1 (0.2-10.0 keV), using PIMMS 
(Mission Count Rate Simulator) Version 3.2, with a power 
law, p hoton index of 2 (as Hasin ger et al. l200ll) . lKuiper et al. I 
(2002a) found 51 sources in a comparable field of view. We 
detected 14 so urces within the cen tral 5' radius (excluding the 
pulsarl wherelMineo et al. I fcOOOh found 7, using the ROSAT 
data of iKuiper et al. I (|1998). The whole field of view can be 
seen in Fig. [4] 



5. Spectral analysis 

The total energy spectrum of the MSP PSR J02 18+4232 can 
be extracted from the MOS imaging data. However, as the 
lightcurves have been seen to vary dramatically from energy 
band to energy band (see Sect. [3}, we have also used the PN 
timing data from which the phase resolved spectral information 
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can be extracted, to isolate the spectra from different regions in 
the lightcurve. We have fitted these spectra with a variety of 
models, in an attempt to find which is the most appropriate. 

5.1. MOS spectral analysis 

We extracted the MOS spectra using an extraction radius of ~1' 
and rebinned the data into 15 eV bins. We used a similar neigh- 
bouring surface, free from X-ray sources to extract a back- 
ground file. We used the SAS tasks 'rmfgen' and 'arfgen' to 
generate a 'redistribution matrix file' and an 'ancillary response 
file' . We then used Xspec (Version 1 1 . 1 .0) to fit the spectra. We 
binned up the data to contain at least 50 counts/bin. Fitting the 
MOS data together, we find that a good fit for the spectrum (be- 
tween 0.2-10.0 keV) is a single absorbed power law, photon in- 
dex of 1 . 10+0.06,*;;= 1 00 (29 degrees o f freedom (dof)) whe n 
the Nj| was frozen at 5xl0 20 cirr 2 , see IVerbunt et al. I (11996). 
Allowing the Njj to vary, gives Njj=(7.6+4.3) x 10 20 crrT 2 
and a photon index of 1.19+0.12,^=1.04 (28 dof). A black- 
body plus a power law also provides a good fit to the data (see 
Table 13 where the unabsorbed flux (ergs cm _2 s _1 ) in the 0.6- 
10.0 keV range can also be found). A single blackbody gives a 
very poor fit to the data, kT=0.85+0.04 keV, ^=3.79 (29 dof), 
when the Njj was frozen at 5xl0 20 cirT 2 . 



5.2. PN spectral analysis 
5.2.1 . The total pulsar spectrum 

The PN spectra were extracted using regions as described in 
Sect. We used the PN 'timing and burst redistribution ma- 
trix file' supplied by the XMM-Newton Science Operations 
Centre 1 . Again we used Xspec to fit the spectra and we binned 
the data to have at least 50 counts/bin. The PN spectrum, be- 
tween 0.6-10.0 keV, is well fitted by a single absorbed power 
law of photon index 1.08 ±0.06 (see Fig.|5j left hand side), in 
good agreement with the MOS data (1.10±0.06). The photon 
index of the power law fit is t he same, within the errors, as that 
found bv lMineo etai~ld2000l) (0.94±0.22), who fitted the spec- 
trum of PSR J02 18+4232, observed by BeppoSAX, between 
2-10 keV. 

The advantage of the PN data is that we can fit the spec- 
tra from different regions in the lightcurve. In all the following 
fits, we assume Njj = 5 xl0 20 crrT 2 and we simply subtracted a 
fraction of the background spectrum (as above) corresponding 
to the fraction of the pulsar phase interval that we were inter- 
ested in, unless stated otherwise. 



5.2.2. The pulsed region 

We first extracted the spectrum of the whole of the pulsed re- 
gion (0.58 < (p < 1.45, see Fig. [2) and binned it up to have at 
least 30 counts/bin. The results of the best spectral fits and the 
unabsorbed flux of this pulse can be found in Table|2] We then 
subtracte d the non-pu lsed region from the pulsed region data, 
as M ineo et al. I .200 0). so that we could study the pulsed emis- 



To+al Source 
(minus BKG) 



1 




10 



Total Pulsed 
(minus BKG & DC) 




10 1 



• 



Energy (keV) 



Energy (keV) 



http://xmm.vilspa.esa.es/external/xmm_sw_cal/calib_frame.shtml 



Fig. 5. The left hand plot shows the power law fit (solid line) to 
the unfolded total (pulsed + non-pulsed) pulsar spectrum using 
EPIC PN data. On the right hand side, the unfolded spectrum 
from the pulsed region only and the best fit to these data are 
shown. The thick dashed line is the best fit to the data above 2 
keV, see Section l5.2.2l and the thin dashed line is the extrap- 
olation of this fit, below 2 keV. Please note that the data have 
been rebinned to contain fewer bins than the data presented in 
Table |2 to make the plots clearer. The best fit is, however, the 
same as that presented in Tabled 



sion only. We find, that a power law gives a good fit to the data. 
The spectrum with the best power law fit, can be seen in Fig. [5] 
right hand side. However, we find that the power law photon in- 
de x of the pulsed em ission (1.19+0.06) is softer than that found 
bv lMineo et ai"1 d200dl) (0.61+0.32), who fitted the BeppoSAX 
data between ~2-10 keV, even taking into account the errors. 
The main difference between the XMM-Newton data and the 
BeppoSAX data is that XMM-Newton has spectral informa- 
tion down to ~0.6 keV, whereas the BeppoSAX spectral infor- 
mation has a lower limit of about 2 ke V. This discrepancy in the 
photon indices of the power law fits to the data sets may there- 
fore stem from an excess of low energy photons in the XMM- 
Newton data, not observable by BeppoSAX. We t herefore tried 
to fit o ur own data in the same energy band as iMineo et al. I 
(2000). Using only this energy band we found that a power 
law, with a photon index of 0.68+0.32 (x]-Q.99, 15 dof), gave 
a good fit to the data. This power law fit can be seen as the thick 
dashed line in Fig. [5] right hand side. The photon index is con- 
sistent with the power law photon index found by M ineo et al. I 
(2000). The thin dashed line indicates the extrapolation of this 
power law fit to lower energies. However, the data points at 
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Table 2. Best fitting models to spectra from the MOS and PN data. In all the fits where no instrument is mentioned, the data is 
from the PN. Regions fitted are the whole lightcurve, the whole of the pulsed region, each of the pulses individually, the interpulse 
and the non-pulsed region (Njj=5 xl0 20 cirr 2 ). The upper part of the table shows the data minus the background. The lower part 
of the table shows the data minus the non-pulsed region and the background. The flux given is unabsorbed flux (ergs cm _2 s _1 ) in 
the 0.6-10.0 keV range, averaged over the whole of the pulse period, with errors of the order +10%. 
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low energies (0.6-2 keV) deviate from the extrapolation of the 
power law fit. Including a blackbody model, where such ther- 
mal emission may be ex pected following theoretical work e.g. 
IZhang &Chengl l20031. gives an equally good fit to the data 
as a single power law fit. We also find that the power law pho- 
ton index is now c ompatible with the photon index found by 
iMineo et al~l d2000l> . with the BeppoSAX data. A two com- 
ponent blackbody model also gives a good fit to the data, see 
TableEl 



If we assume that the spectru m above 2 k eV is best fitted by 
the power law given in lMineo et al. 1 ll2000l) (as we have ruled 
out the single blackbody fit and also shown that we find a sim- 
ilar fit for our data above 2 keV), we can use this value to try 
to determine whether the single power law or a blackbody plus 
a power law gives the better fit to the data . Fitting the pulsed 
spectrum with the power law fixed at the Min eo et al. I (|2000) 
value, we find ;^;=1.79 (43 dof) and doing the same for the 
two component model, we find ^=1.17 (41 dof) (blackbody, 
kT=0. 18+0.04). Using an F-test we found that the improve- 
ment in the fit was significant at the 3.8<x level. Even if we 
used the softes t power law s l ope p ossible, taking into account 
the error on the Mi neo et al. 1 ( 2000) data, we found that the im- 
provement in the fit was significant at the 2.6<x level. Thus we 
may have some indication that the two component model may 
be the best fit to the data. 



5.2.3. Pulse 1 



We then extracted the spectrum from pulse 1 only (0.1 1 < <p < 
0.45), binned it up to have at least 30 counts/bin and subtracted 
the background. The results of the best spectral fits and the un- 
absorbed flux of this pulse can be found in Table [2] Using the 
F-test we find that all the models that we have tried are equally 
probable as the most appropriate fit. The data from pulse 1 
plotted with the best power law fit can be found in Fig.[6](left 
hand side). A single blackbody and a thermal bremsstrahlung 
model give poor fits to the data. Subtracting the non-pulsed re- 
gion from the pulse 1 data, as before, we find similar spectral 
fits to those found with the simple background extraction. The 
data plotted with the best power law fit can be seen in Fig. [6] 
right hand side. The power law spectral index that we find us- 
ing the XMM-Newton data b etween 0.6-10 keV (1.15+0.13) 
is larger than that found bv lMineo et"aD d2000h (2-10 keV), 
who found 0.84+0.35, although it is compatible within the er- 
ror bars. We did, however, try to fit our data above 2 keV, as in 
Section 15.2.21 We found that a harder power law with photon 
index 0.66+0.29 gave a good fit to the data, X 2 ..=\-29 (1 9 dof). 
This is again compatible with the Mi neo et al. I ( 120001) result. 
The power law fit above 2 keV can be seen as the thick dashed 
line in Fig.|6] right hand side. The thin dashed line indicates the 
extrapolation of this power law fit to lower energies. 
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Fig. 6. Plots showing the power law fit (solid line) to pulse 1. 
On the left is pulse 1 minus the background contribution and on 
the right is pulse 1 minus the DC and background contributions. 
The thick dashed line is the best fit to the data above 2 keV, see 
Section l5.2.3l and the thin dashed line is the extrapolation of 
this fit to below 2 keV. These data are not averaged over the 
whole of the period, so the exposure time = Ac/). T (where A<p is 
the phase over which the pulse was extracted). Please note that 
the data have been rebinned to contain fewer bins than the data 
presented in Tabled to make the plots clearer. The best fit is, 
however, the same as that presented in Tabled 

5.2.4. Pulse 2 

Extracting the spectrum from pulse 2 only (0.58 < <p < 0.88), 
binning it up to have at least 20 counts/bin and subtracting the 
background, we find the fits presented in Table |2 Again a sin- 
gle blackbody and a thermal bremsstrahlung model give poor 
fits to the data. The data plotted with the power law model can 
be found in Fig. (left hand side). We subtracted the non- 
pulsed data from the pulse 2 data (as above). We again find 
similar results for the power law fit, as with the simple back- 
ground extraction, however, as we found with the pulsed emis- 
sion, the power law photon index is softer than that found by 
iMineo et al. I ((2000) (0.42+0.36), again taking into account the 
errors. These data and the power law model can be found in 
Fig.Q(right hand side). 

5.2.5. The interpulse region 

Fitting the spectrum to the interpulse region (0.88 < (p < 1.11) 
in a similar way, we find the fits as given in Table[2] The errors 
are much larger in this case, due to the small number of bins. 



Fig. 7. Plots showing the best fit (power law, solid line) to pulse 
2. On the left is pulse 2 minus the background contribution and 
on the right is pulse 2 minus the DC and background contribu- 
tions. These data are not averaged over the whole of the period, 
so the exposure time = A<p. T (where A<f> is the phase over which 
the pulse was extracted). Please note that the data have been re- 
binned to contain fewer bins than the data presented in Table|2] 
to make the plots clearer. The best fit is, however, the same as 
that presented in Tabled 

Again a single blackbody and a thermal bremsstrahlung model 
give poor fits to the data. We subtracted the non-pulsed data 
from the interpulse data (as above), however this resulted in 
too few counts to create a spectrum. 

5.2.6. The non-pulsed region 

Fitting the spectrum of the non-pulsed region (0.45 < <p < 0.58) 
in a similar way (results in Table|2j, we find that a simple power 
law gives the best fit to the data. In this case, a blackbody plus 
a power law model does not give a good fit. Taking the non- 
pulsed flux (0.6-10.0 keV) given in Table[2] and the total flux 
(0.6-10.0 keV) also given in Table|5] we find a pulsed fraction 
of 61*5%, consistent with the pulsed fraction given in Tabled 
(0.6-12.0 keV). 

6. Discussion 

Than ks to the la rge collecting area of XMM-Newton 
djansen et al. I l200ll) . even from a short observation of PSR 
J0218+4232, we are able to undertake detailed spectral anal- 
ysis of the various components of the lightcurve. 
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PSR J0218+4232 shows a very hard spectrum and an 
analogy is often made between this MSP and two others, 
PSR B 182 1-24 and PSR B 1937+21, which also show pulsed 
emission, large E values (2-20 x 10 35 er gs s~') and X-ray 
spectra with photon i ndices less than 2 dSaito et al~l 1 19971 
lTakahashietaLll200lh . These pulsars are in turn compared to 
the Crab pulsa r, which shows narrow pulsations and a hard X- 
ray spectrum jToor & Seward ill 9771) . It is believed that such 
hard emission arises from the magnetosphere. However, the 
Crab pulsar shows pulses which have FWHM of ~0.042 and 
0.075 in phase at similar X-ray energies dEikenberrv & Faziol 
119971) . which are 2-3 times narrower than those of PSR 
J02 18+4232. PSR B 1 82 1 -24 also has a pulse widt h of less than 
0.04 in phase jRots et al. ll998HSaito et al. 1 19971) and those of 
PSR B 1937+21 are also only -0.06 in phase (|Takahashi et al. I 
l200ll) . The larger width of the PSR J02 18+42 32 pulses may 
simpl y be due t o it being an ali gned rotator jNavarro et al. I 
Il995l) . However, iLuo et al I J20001) state that with their polar 
cap model, they can reproduce the wide double-peaked profile 
without assuming an aligned rotator. No netheless, the particu- 
larly hard spectrum, extending to y-rays dKuiper et al. 1 12002b. 
2000) and the wider pulses shown by PSR J0218+4232 could 
indicate that it is a different kind of MSP to PSR B 182 1-24 and 
PSRB1937+21. 



6.1. The pulsed emission 



We have seen that the form of the pulsed emission varies sig- 
nificantly with energy (see Sect.[3]and Fig.[3J. We have there- 
fore tried a variety of simple models to fit the pulsed emission 
data, see Sect.|5]and Table|2] For the pulsed spectrum and the 
spectra of pulses 1 and 2, we find that a single blackbody fit 
is not a good description of the data, but that a single power 
law does give a good fit. However, the photon indices that 
we find when fitting the XMM- Newton data, are significantly 
softer than those found by Mineo et al. when fitting the 

spectra observed with the BeppoSAX MECS between 2 and 
10 keV. By fitting the XMM-Newton spectra with a two com- 
ponent model (black body plus power la w, in the same way as 
for PSR J0437-4715 JZavlin et al. 120021)') we recover th e same 
power law photon indices found by M ineo et al. (2000). which 
could indicate that there may be additional emission below 2 
keV and therefore not observable with BeppoSAX. Taking the 
blackbody temperature (2.90+0.70 xlO 5 K) from the black- 
body plus power law fit to pulse 1, we find a radius of emission 
of 0.37±0.33 km (90% confidence limit), consistent with the 
tem perature and radius of a polar cap (10 -10 7 K a nd ~1 km 
e.g. lzhang & Cheng 1120031: IZavlin & Pavlov 111 998l and refer- 
ences therein). 

There is some evidence from both the spectral fitting and 
the lightcurves presented in Sect. [5] that the interpulse region 
emit s a softer spectrum than that of the two pulses (as also seen 
by Kuine r et al. 120 02a). 



6.2. The non-pulsed emission 

With regards to the non-pulsed region, we find that a sim- 
ple blackbody, with a temperature consistent with that deter- 
mined using t he ROSAT HRI, Bep poSAX MECS and Chandra 
instruments (Kuip eret al. l E002a). can fit the data. However, 
the radius of emission determined from this blackbody is only 
0.10+0.08 km. This indicates that from the small radiu s, the 
emissio n is unlikely to be from a polar cap. iKuiper et al. I 
(2002a) suggest that this could nonetheless be thermal emission 
from the polar caps, if PSR J0218+4232 is a nearly aligned ro- 
tator and the viewing angle is small, so that the polar caps are 
always in sight. Alternatively, the data can be well described 
by a single power law, where the photon index is the same, 
within the errors, as is found to fit all the other sections of 
the lightcurve. This would be consistent with re cent polar cap 
models, such as that of Zhang & Cheng N2003I) . where X-ray 
emission in these models consists of one power law and two 
thermal components. In these cases the nonthermal X-rays are 
produced by the synchrotron radiation of e* pairs created in 
a strong magnetic field near the stellar surface; the soft ther- 
mal X-rays are produced by heating of polar cap areas with a 
radius of the order of 1 km; and the medium hard thermal X- 
rays result from the polar cap heating by the return particles 

from the outer gap. This is consistent with data taken at higher 

I n r 1 1 1 1 1 I I 

energies (see Kuip er et al. 1 12000. 2002b), where Kuiper et al. I 

(2000) determine that extrapolating a hard power law spectrum 

from the XMM-Newton spectral range is just in agreement with 

the OSSE upper limit. 

7. Conclusions 

We have presented XMM-Newton MOS imaging (0.2-10.0 
keV) and PN timing (0.4-12.0 keV) data of PSR J0218+4232. 
Using the large number of counts, we have confirmed the 
previously detected pulsations of PSR J0218+4232 and we 
have shown that the form of the lightcurve varies with energy. 
We have shown that the broad band X-ray spectrum between 
0.6-10.0 keV of the pulsed emission is softer than the spec- 
trum measured between 2-10 keV by BeppoSAX, and we have 
shown some evidence that a simple power law may not be a 
good description of the 0.6-10.0 keV spectrum. A longer ob- 
servation of this millisecond pulsar should establish whether or 
not this is indeed the case. 
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